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(57) ABSTRACT 

A method, radio receiver, and system to autonomously 
receive and decode a plurality of signals having a variety of 
signal types without a priori knowledge of the defining char- 
acteristics of the signals is disclosed. The radio receiver is 
capable of receiving a signal of an unknown signal type and, 
by estimating one or more defining characteristics of the 
signal, determine the type of signal. The estimated defining 
characteristic(s) is/are utilized to enable the receiver to deter- 
mine other defining characteristics. This in turn, enables the 
receiver, through multiple iterations, to make a maximum- 
likelihood (ML) estimate for each of the defining character- 
istics. After the type of signal is determined by its defining 
characteristics, the receiver selects an appropriate decoder 
from a plurality of decoders to decode the signal. 
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SELF-CONFIGURABLE RADIO RECEIVER 
SYSTEM AND METHOD FOR USE WITH 
SIGNALS WITHOUT PRIOR KNOWLEDGE 
OF SIGNAL DEFINING CHARACTERISTICS 

5 

CROSS-REFERENCE TO RELATED 
APPLICATIONS 

The present application derives priority from U.S. Provi- 
sional Application No. 60/736,880 filed 12 Nov. 2005. 

STATEMENT OF GOVERNMENT INTEREST 

The invention described hereunder was made in the per- 
formance of work under a NASA contract, and is subject to 
the provisions of Public Law #96-517 (35 U.S.C. 202) in 
which the Contractor has elected not to retain title. 

BACKGROUND 20 

a. Field of Invention 

The invention relates to radio receivers and, more particu- 
larly, to an autonomous radio receiver capable of receiving a 
signal without a priori knowledge about the defining charac- 25 
teristics of the signal, recognizing these characteristics 
directly from the received signal, and appropriately decoding 
the signal based on the signal type. 

b. Background of the Invention 

Conventional radio receivers typically have complete prior 30 
knowledge of the defining characteristics of signals they are 
designed to receive. Because conventional radio receivers 
know what type of signal they are going to receive they are 
typically hardwired or manually configured to do so. Con- 
versely, conventional radio receivers are incapable of receiv- 35 
ing signal types different from that which they were intended 
for. In some cases where they do have such capability, con- 
ventional radio receivers require specific pre-configuration 
according to a pre-determined schedule. Thus, conventional 
radio receivers are capable of only receiving pre-deteimined/ 40 
pre-programmed types of signals. For example, cellular tele- 
phones are designed according to several different standards, 
and single-mode telephones are typically preconfigured to 
receive only one type of cellular telephony signal. Indeed, a 
single-mode telephone configured to use Code Division Mul - 45 
tiple Access (a popular system in the United States) is unable 
to receive a signal using Global Systems Mobile (GSM) sig- 
nals, which is the European standard. While there are dual- 
mode telephones to solve the CDMA/GSM dichotomy, there 
remains the issue and potential future issue of varying fre- 50 
quency bands, data rates, modulation types, modulation indi- 
ces, pulse shapes, and so on, which may affect the usability of 
cellular telephones. Therefore, there is a need for a system 
capable of receiving multiple signal types without a priori 
knowledge about the defining characteristics of the signals, 55 
and accurately recognizing defining features of the incoming 
signals and respond intelligently without explicit pre-con- 
figuration or reprogramming to define the functions of the 
receiver. In other words, there is a need for a system capable 
of recognizing multiple signal types and autonomously 60 
reconfiguring itself to decode the particular type of signal that 
is received. 

SUMMARY OF THE INVENTION 

65 

It is, therefore, an object of the present invention to provide 
a system for autonomously receiving and decoding a plurality 
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of signals having a variety of signal types without prior 
knowledge of the defining characteristics of the signals. 

According to the present invention, the above-described 
and other objects are accomplished by providing a system 
inclusive of hardware and a suite of software modules to 
autonomously recognize various defining characteristics of a 
signal that are unknown prior to reception. Some of these 
defining characteristics include, but are not limited to, the 
angle of arrival, data rate, symbol timing, carrier frequency 
and phase, modulation index, modulation type, signal-to- 
noise ratio (SNR), code type, and decoded message bits. The 
signal characteristics are recognized by estimating one or 
more characteristics of the signal, and the estimated charac- 
teristic is utilized to enable the apparatus and system to deter- 
mine other defining characteristics. Through multiple itera- 
tions of estimation and calculation, the system is able to make 
a maximum-likelihood (ML) or near-ML estimate for each of 
the defining characteristics. The ML estimates assist in iden- 
tifying the type of signal and the system selects an appropriate 
decoder from a plurality of decoders to eventually decode the 
signal. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Additional aspects of the present invention will become 
evident upon reviewing the embodiments described in the 
specification and the claims taken in conjunction with the 
accompanying figures, wherein like numerals designate like 
elements, and wherein: 

FIG. 1 is a block diagram of a radio receiver 100 operating 
in a coarse mode for determining the signal type of an incom- 
ing signal with unknown defining characteristics in accor- 
dance with one exemplary embodiment of the invention. 

FIG. 2 is a block diagram of radio receiver 100 operating in 
a fine mode in accordance with one exemplary embodiment 
of the invention. 

FIG. 3 is a block diagram of the detailed operations of radio 
receiver 100. 

FIG. 4 is a block diagram of a system 200 for receiving and 
decoding a signal without a priori knowledge of the defining 
characteristics of an incoming signal. 

FIG. 5 is a block diagram of an exemplary embodiment of 
modulation index estimator 1110. 

FIG. 6 is a block diagram of one exemplary embodiment of 
closed-loop frequency estimator 1220. 

FIG. 7 is a block diagram illustrating one exemplary 
embodiment of a frequency correction error signal detector 
1225' for a residual carrier signal. 

FIG. 8 is a block diagram of an exemplary embodiment of 
a frequency correction error signal detector 1225" for a 
residual carrier signal. 

FIG. 9 is a block diagram of a closed-loop frequency esti- 
mator 1230 for a residual carrier signal, with the approximate 
error signal detector. 

FIG. 10 is a block diagram of a circuit for a digital loop 
filter 1250 for frequency tracking loops. 

FIG. 11 is a block diagram of an open-loop frequency 
estimator 1260 for a suppressed carrier signal with known 
data. 

FIG. 12 is a block diagram of an error signal detector 1275 
for the closed-loop frequency correction for a suppressed 
carrier signal with known data. 

FIG. 13 is a block diagram of one embodiment of a closed- 
loop frequency estimator 1270 for a suppressed carrier signal 
with known data. 
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FIG. 14 is a block diagram of a closed-loop frequency 
estimator 1280 for a suppressed carrier signal with known 
data. 

FIG. 15 is an all-digital version of a closed-loop frequency 
estimator 1290 for a suppressed carrier signal with known 
data. 

FIG. 16 is a block diagram of an open-loop frequency 
estimator 1292 for a suppressed carrier signal with unknown 
data. 

FIG. 17 is a block diagram of a frequency correction error 
signal detector 1294 for a suppressed carrier signal with 
unknown data. 

FIG. 18 is a block diagram of a closed-loop frequency 
estimator 1295 for a suppressed carrier signal with unknown 
data. 

FIG. 19 is a block diagram of a closed-loop frequency 
estimator 1296 for a suppressed carrier signal with unknown 
data. 

FIG. 20 is a block diagram of one embodiment of a low 
SNR data format classifier 1321 for a binary phase-shift key- 
ing (BPSK) signal. 

FIG. 21 is a block diagram of an embodiment of a high 
SNR data format classifier 1322 for a BPSK signal. 

FIG. 22 is a block diagram of a low SNR coherent data 
format classifier 1323 for quadrature phase-shift keying 
(QPSK) modulation. 

FIG. 23 is a block diagram of a reduced complexity low 
SNR data format classifier 1324 for BPSK modulation. 

FIG. 24 is a block diagram of one exemplary embodiment 
of a reduced complexity high SNR data format classifier 1325 
for BPSK modulation. 

FIG. 25 is a block diagram illustrating an embodiment of 
SSME SNR estimator 1341 for M-ary phase shift keying 
(MPSK) modulation, sampled version. 

FIG. 26 is a block diagram of a complex baseband SSME 
SNR estimator 1342 for MPSK modulation, integrate and 
dump version. 

FIG. 27 is a tree diagram for the sampled version of the 
SSME SNR estimator 1341 for MPSK modulation. 

FIG. 28 is a block diagram of one exemplary embodiment 
of time-multiplexed SSME SNR estimator 1346 for MPSK 
modulation. 

FIG. 29 is a flow diagram of one exemplary embodiment of 
an adaptive scheme 1344 for SNR estimation of MPSK 
modulation. 

FIG. 30 is a block diagram of a system to finely integrate 
and sample a continuous time signal to obtain a high rate 
discrete-time signal. 

FIG. 31 is a digital implementation of a SSME-based data 
rate estimation system 1362. 

FIG. 32 is a block diagram of one exemplary embodiment 
of an MPSK modulation classifier 1410. 

FIG. 33 is a block diagram of a high SNR data transition 
tracking loop (DTTL) 1620. 

FIG. 34 is a block diagram of a MAP-motivated linear 
DTTL 1630. 

FIG. 35 is a block diagram of a maximum likelihood slid- 
ing-window estimator of symbol timing for an NRZ data 
stream. 

FIG. 36 is a block diagram of a weighted running accumu- 
lator. 

FIG. 37 is a block diagram of a weighted running accumu- 
lator 1650. 
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FIG. 38 is block diagram of a real non-coherent version of 
a cross -spectrum synchronizer 1660. 

FIG. 39 is a block diagram illustrating NC-LDTTL 1670. 

5 DETAILED DESCRIPTION OF THE PREFERRED 
EMBODIMENT 

The present invention is an autonomous radio receiver 
capable of receiving a signal without prior knowledge about 
to the defining characteristics of the signal, instead recognizing 
these characteristics directly from the received signal, and 
appropriately decoding the signal based on the signal type. 

FIG. 1 is a block diagram of a radio receiver 100 for 
determining the signal type of an incoming signal with 
15 unknown defining characteristics in accordance with one 
exemplary embodiment of the invention. Radio receiver 100 
includes a modulation index module (MIM) 110 that may be 
any hardware and/or software suitably configured to estimate 
the modulation index and/or determine the actual modulation 
20 index of an incoming signal with unknown characteristics. It 
is common is to allocate a portion of the total transmitted 
power to a discrete carrier for purposes of carrier synchro- 
nization. For example, in the case of binary phase-shift key- 
ing (BPSK) modulation, the simplest way to accomplish this 
25 is to employ a phase modulation index other than 90 degrees . 
When this is done, the fraction of power allocated to the 
discrete carrier P c becomes a function of the total transmitted 
power P r MIM 110 estimates and/or determines the fraction 
of total power that is allocated to a non-modulated carrier 
30 signal and to a data -modulated signal and outputs an esti- 
mated modulation index including such estimate. Total power 
is the sum of the non-modulated carrier signal power and the 
data-modulated signal power, in watts. 

The radio receiver 100 also includes a frequency correction 
35 module (FCM) 120 in communication with MIM 110. FCM 
120 may be any hardware and/or software suitably configured 
to track the offset frequency of the incoming signal with 
unknown characteristics without knowledge a priori of the 
signal-to -noise ratio (SNR), data rate, modulation type, etc. 
40 Thus, FCM 120 receives the incoming signal and the esti- 
mated modulation index from MIM 110, and outputs a signal 
having an estimated offset frequency based on the incoming 
signal and the estimated modulation index. FCM 120 is in 
communication with a multiplier 125 and outputs the signal 
45 having an estimated frequency offset to multiplier 125, and 
multiplier 125 mixes the original signal with the signal having 
the estimated frequency offset to produce a frequency-cor- 
rected received signal. 

Radio receiver 100 also includes a joint estimator module 
50 (JEM) 130 in communication with MIM 110 and multiplier 
125. JEM 130 may be any hardware and/or software suitably 
configured to receive the estimate of the modulation index 
and/or the actual modulation index from MIM 110 and to 
utilize the estimated modulation index and frequency cor- 
55 rected received signal to estimate and/or determine defining 
characteristic(s) of the incoming signal. In the preferred 
embodiment JEM 130 utilizes the estimated modulation 
index to estimate and/or determine the data rate of the incom- 
ing signal, the SNR of the incoming signal, the waveform 
60 shape and pulse shape of the incoming signal, and/or symbol 
timing. 

As discussed in greater detail below, one or more of the 
estimates and/or characteristics produced by JEM 130 is 
transmitted back to MIM 110 to refine the estimate of the 
65 modulation index and/or help determine the modulation 
index. Furthermore, one or more of the estimates and/or char- 
acteristics produced by JEM 130 is transmitted to a modula- 
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tion classifier module (MCM) 140, a carrier synchronization 
module (CSM) 150, and/or a fine symbol synchronizer (FSS) 
160, all in communication with JEM 130 and with each other. 

MCM 140 may be any hardware and/or software suitably 
configured to derive approximations of the maximum likeli- 
hood (ML) modulation classification based on one or more 
estimations and/or characteristics received from multiplier 
125 (i.e., MIM 110 and FCM 120 estimates), JEM 130, CSM 
150, and/or FSS 160. 

CSM 150 may be any hardware and/or software suitably 
configured to track the carrier phase independent of the 
modulation of the incoming signal. Furthermore, CSM 150 is 
such that it includes a structure that performs the carrier 
synchronization function for all orders of phase-shift-keying 
(PSK) modulation based on one or more estimations and/or 
characteristics received from MIM 110, multiplier 125 (i.e., 
MIM 110 and FCM 120 estimates), JEM 130, MCM 140, 
and/or FSS 160. 

FSS 160 may be any hardware and/or software suitably 
configured to determine an estimate of the symbol timing of 
the incoming signal. FSS 160 determines the estimate of the 
symbol timing of the incoming signal based on one or more 
estimates and/or characteristics received from signal com- 
biner 125 (i.e., MIM 110 and FCM 120 estimates), JEM 130, 
and/or CSM 150. 

In general operation, radio receiver 100 functions to 
receive an incoming signal, and through estimations and one 
or more refining iterations, determines the signal type and 
decodes the incoming signal utilizing an appropriate decoder. 
How MIM 110, FCM 120, JEM 130, MCM 140, CSM 150, 
and FSS 160 calculate the estimations and perform the itera- 
tions to determine the signal type is discussed in greater detail 
below. Notably, the workings of MIM 110, FCM 120, JEM 
130, MCM 140, CSM 150, and FSS 160 depend on whether 
radio receiver 100 is performing coarse estimations, fine esti- 
mations, non-coherent estimations, or coherent estimations, 
and on what information has been estimated and/or deter- 
mined by prior iterations. 

Again, MIM 110 is configured to receive signals with 
unknown characteristics. Because radio receiver 100 does not 
know exactly what type of signal it may receive, MIM 110 is 
configured to receive and estimate the modulation index for 
an array of signal types. 

For example, phase-shift keying (PSK) is a method of 
digital communication in which the phase of a transmitted 
signal is varied to convey information. There are several 
methods that can be used to accomplish PSK. The simplest 
PSK technique is called binary phase-shift keying (BPSK). It 
uses two opposite signal phases (0 and 180 degrees). In the 
case of binary phase-shift keying (BPSK) modulation, a 
phase modulation index (3 other than 90 degrees may be used. 
When this is done, the fraction of power allocated to the 
discrete carrier becomes P C =P, cos 2 |3 with the remaining 
fractional power P^=P, sin 2 p available for data detection. 
When using this signaling mode, the power spectrum of the 
data modulation should be such that it does not interfere with 
the extraction of the discrete carrier by an appropriate carrier 
tracking loop such as a phase-locked loop (PLL). In other 
words, the discrete carrier should be inserted at a point where 
the power spectrum of the data modulation is minimum, and 
preferably equal to zero. In the case of digital data, this rules 
out direct modulation of the carrier with a Non-Return to Zero 
(NRZ) data stream whose spectrum is maximum at direct 
current (DC) which at radio frequencies would correspond to 
the carrier frequency. Thus, for example, the data may be first 
modulated onto a sub-carrier whose frequency is selected 
significantly higher than the data rate so that the sidebands of 
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the data modulation are sufficiently reduced by the time they 
reach the carrier frequency. More sophisticated forms of PSK 
exist. In multiple phase-shift keying with or without phase 
compensation. In the most general scenario, a taxonomy of 
5 cases for analysis is illustrated by the tree diagram of FIG. 27 
for the sampled version of the SSME 1341 . FIG. 27 begins at 
the square node in the middle and proceeds outward to any of 
the eight leaf nodes representing combinations of co, oo, co^, 
and ct^. For SSME 1342, a few of the tree branches of FIG. 27, 
namely, 2c and 3c, do not apply. 

Other signals use a coherent communication mode where 
carrier synchronization is established directly from the data- 
bearing signal (e.g., using a conventional Costas loop). In this 
15 case, none of the transmitted power is allocated to a discrete 
carrier and the system is said to operate in a suppressed carrier 
mode which, in the case of BPSK, corresponds to (3=90°. 
Although a Costas loop operates with a less efficient perfor- 
mance (e.g., larger mean-square phase tracking error) than a 
20 PLL, it offers the advantage of not requiring a discrete carrier 
to lock onto and thus all of the transmitted power can be used 
for the purpose of data detection. 

Given that the foregoing or other transmission modes may 
be received by radio receiver 100, it is important that MIM 
25 1 10 be able to estimate the modulation index or equivalently, 
the ratio of transmitted carrier to data power of the various 
received signals. In one exemplary embodiment, MIM 110 
utilizes an ML estimation approach to estimate the modula- 
tion index along with appropriate approximation of the non- 
30 linearities that result, to allow for practical implementations 
at course and fine modes. 

With reference again to FIG. 1, radio receiver 100 is shown 
operating in a coarse mode according to one exemplary 
embodiment of the invention. When radio receiver 100 
35 receives an incoming signal (which is represented as r(t)), the 
signal is transmitted to MIM 110 and signal combiner 125. 
MIM 100 estimates an initial modulation index without a 
priori knowledge of the modulation, data rate, and/or symbol 
timing of the incoming signal. The incoming signal with the 
40 initial estimate of the modulation index is then transmitted to 
other modules of radio receiver 1 00 to estimate other defining 
characteristics. 

If the signal with the initially estimated modulation index 
includes a residual carrier, the signal is transmitted to FCM 
45 120 for frequency acquisition and tracking. If the signal with 
the initial estimated modulation index does not include a 
residual carrier, the signal is also transmitted to FCM 120 for 
frequency acquisition and tracking. Furthermore, the signal 
with the initially estimated modulation index and corrected 
50 frequency is also transmitted to JEM 130, 140, CSM 150 
AND/OR 160. 

The signal with the initially estimated modulation index is 
transmitted to JEM 13 0 so that the received signal can be used 
to estimate an initial data format (via data format module 
55 132), an initial SNR (via SNR module 134), an initial data rate 
(via data rate module 136), and/or an initial estimate of the 
symbol timing (via CSS 138). Furthermore, the signal with 
the initial estimated modulation index is transmitted to CSM 
1 50 so that the received signal can be used to initially track the 
60 carrier phase of the incoming signal. Moreover, CSM 150 
also receives and uses the initial estimates of the data format 
(from data format module 132), the SNR (from SNR module 
134), the data rate (from data rate module 136), and/or the 
estimate of the symbol timing (from CSS 138); and/or 
65 receives the original signal combined with the output of F CM 
120 from multiplier 125 in initially tracking the carrier phase 
of the incoming signal. 
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FSS 160 receives and uses the initial estimates of the data 
format (from data format module 132), the SNR (from SNR 
module 134), and/or the data rate (from data rate module 
136); and/or receives the original signal combined with the 
output of FCM 120 from multiplier 125 to initially estimate 5 
the symbol timing of the incoming signal. Likewise, MCM 
140 receives and uses the initial estimates of the data format 
(from data format module 132), the SNR (from SNR module 
134), the data rate (from data rate module 136), and/or the 
estimate of the symbol timing (from CSS 138); and/or 10 
receives the original signal combined with the output of FCM 
120 multiplier 125 to initially estimate the modulation clas- 
sification of the incoming signal. In one exemplary embodi- 
ment, the initial estimate of the carrier phase tracking, the 
initial estimate of the symbol timing, the initial estimate of the 1 5 
modulation classification, the initial estimate of the data for- 
mat, and/or the initial estimate of the data rate are used to 
refine and/or determine the type of the incoming signal in a 
fine mode of operation of radio receiver 100. 

FIG. 2 is a block diagram of radio receiver 100 operating in 20 
a fine mode in accordance with one exemplary embodiment 
of the invention. When operating in fine mode, MIM 110 
receives the initial estimate of the carrier phase tracking (from 
CSM 150), the initial estimate of the symbol timing (from 
FSS 160), the initial estimate of the modulation classification 25 
(from MCM 140), the initial estimate of the data format (from 
data format module 132), and/or the initial estimate of the 
data rate (from data rate module 136) to determine a better 
estimate of the modulation index. This better estimate of the 
modulation index is then utilized similar to the initial estimate 30 
of the modulation index discussed above with respect to FIG. 

1. Thus, at the end of operating in fine mode, radio receiver 
100 will produce a better estimate of the carrier phase track- 
ing, the symbol timing, the modulation classification, the data 
format, and/ or the data rate, in which one or more may be used 35 
in one or more subsequent iterations to produce even finer 
estimations. 

For reference purposes, FIG. 3 is a block diagram of the 
detailed operations of radio receiver 100. The MIM 110 esti- 
mates and/or determines the action of total power that is 40 
allocated to a non-modulated carrier signal and to a data- 
modulated signal and outputs an estimated modulation index 
including such estimate. FCM 120 receives the incoming 
signal and the estimated modulation index from MIM 110, 
and outputs a signal having an estimated corrected frequency 45 
based on the incoming signal and the estimated modulation 
index. FCM 120 outputs the signal having an estimated cor- 
rected frequency to multiplier 125 which mixes the original 
signal with the signal having the estimated frequency correc- 
tion to produce a frequency corrected received signal. JEM 50 
130 receives the modulation index estimate and frequency 
corrected received signal to estimate various defining char- 
acteristic^) of the incoming signal, including data rate of the 
incoming signal, SNR of the incoming signal, waveform 
shape and pulse shape of the incoming signal, and/or symbol 55 
timing. One or more of the estimates and/or characteristics 
produced by JEM 130 is transmitted to MCM 140, CSM 150, 
and FSS 160 (CSM 150 to track the carrier phase, FSS 160 to 
estimate of the symbol timing of the incoming signal, and 
MCM 140 for maximum likelihood (ML) modulation classi- 60 
fication or the modulation classifier. Radio receiver 100 func- 
tions to receive an incoming signal, and through estimations 
and one or more refining iterations, determines the signal type 
and decoding the incoming signal utilizing an appropriate 
decoder. 65 

FIG. 4 is a block diagram of a system 200 for receiving and 
decoding a signal without a priori knowledge of the defining 
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characteristics of an incoming signal. In one exemplary 
embodiment, system 200 includes a radio receiver 210 simi- 
lar to embodiments of radio receiver 100 discussed above 
connected to and/or in communication with an antenna 220, 
wherein antenna 220 may be any hardware and/or software 
suitably configured to receive communication signals. 

In addition, system 200 includes a decoder module 230 
connected to and/or in communication with radio receiver 
210. In one exemplary embodiment, decoder module 230 
includes a module 234 including hardware and/or software 
suitably configured to receive an output signal (including fine 
estimations of the carrier, phase tracking, the symbol timing, 
the modulation classification, the data format, and/or the data 
rate of an incoming signal) from radio receiver 210 and deter- 
mine the type of signal based on a comparison of the estimates 
with the characteristics of an array of known signal types. 
Furthermore, decoder module 230 includes a plurality of 
decoders 238 to decode the incoming signals based upon the 
determination made by module 234, wherein each decoder of 
the plurality of decoders 238 may include any hardware and / 
or software suitably configured to decode a particular signal 
type. As such, each decoder may be any decoder known in the 
art or developed in the future. 

System 200, in various exemplary embodiments, includes 
a display 240 to display the data decoded by decoder module 
230, wherein display 240 may be any hardware and/or soft- 
ware suitably configured to display the data in the incoming 
signal such that the data can be detected by at least one the five 
senses of a user and/or an electronic device (e.g., a computing 
device). For example, display 240 may be an audio and/or 
visual device such as a television, a monitor, a speaker, a 
radio, a computer, and the like. 

The following are examples of the estimation algorithms 
incorporated in the above-referenced modules. 

Example 1 

MIM 110 Coherent Estimation of Modulation Index 
For Binary Phase-Shift Keying (PSK) Signals 

A received BPSK signal r(t) may be represented by the 
following equation (1): 


r(r) = V2P f sinl co c t + /? ^ c n p(t - nT ) + n(t ) 


= V2 P t cos/?sin oj c t + yj2P t sin ft c n p(t - nT)cosco c t + n(t ) 


= V2pT sinw c r + yj2P d ^ c n p(t - nT)cosa) c t + n(t ) 


where p(t) is a pulse shape taking on values ±1 and satis- 
fying T _1 / 0 r p 2 (t)dt=l, T is the symbol rate, {c„} is a binary 
sequence which is treated as independent, identically distrib- 
uted (IID) data taking on values ±1 with equal probability, (I 
is the modulation index, oo c is the carrier frequency in rad/s, 
and n(t) is a bandpass additive white Gaussian noise (AWGN) 
source with single-sided power spectral density N 0 W/Hz. 
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Based on the above AWGN model, then for an observation 
of K data intervals, the maximum likelihood (ML) modula- 
tion index estimation may be derived. Suitable approxima- 
tions to the nonlinearity can be made, corresponding to low 
and high data detection SNR conditions . For high data detec - 5 
tion SNR, the ML modulation index estimation after approxi- 
mation is given by 


cot p - 



r(t)sin(x) c t d t 


£ K +l)T \r(t)p(t-kT)cos(x) c tdt\ 
k = o KI 


( 2 ) 10 


Applying the approximations in (8) to (6) results in the ML 
estimate: 


of 


N 0 I r(t)sma) c tdt 


V" 2 M/2-1 „ 

> 77 2 4(9) 
M q =0 

k = 0 


(9) 


Example 2 


For low data detection SNR, the ML modulation index 
estimation after approximation is given by 


MIM Non-Coherent Estimation of Modulation Index 
For PSK Signals 


cot ft = - 


N 0 f r(t)sm(x) c tdt 


(3) 


20 


2V2U 2 j[ , ‘ +1)7 |r(f)p(r - kT)cosio c tdt\ 


Note that in this case a solution to this equation may not 25 
always exist. 

For M-PSK modulation (M>2), the received signal r(t) can 
be represented in passband by: 


In another embodiment, MIM 110 performs modulation 
index estimation for the case where carrier synchronization 
has not yet been established (i.e., where the carrier phase is 
random). In the non-coherent case, the modulation index 
estimate is formed in the absence of carrier synchronization. 
For BPSK modulation, the received signal is again modeled 
as in equation (1) with the addition of an unknown (but 
assumed to be uniformly distributed) carrier phase to both the 
discrete and data-modulated carriers. This allows a non-co- 
herent estimator of modulation index: 


30 


r(t) = ^f2PfsmoJ c t+^f2Pf cosft 


where 0 w =[2q 77 +(l+(-l) 7Wf/2 )/2] Jt/M is the data modulation 35 
for the nth MPSK symbol, with independent and uniformly 
distributed q, ? . {0, 1, . . . , M-l}. Differentiating the condi- 
tional likelihood of r(t) given 0 W and |j with respect to (3 and 
equating the result to zero results in the transcendental equa- 
tion: 40 


f 


(M/2)-l 

k-i Z 


r(t)sin(i) c tdt = (cot 




2 V 2 P t sin/? 

“/vT 


Xk(q) 


(6) 


45 


*=o j2y/2Pfsmp 


V J 2 

> coshl - 


No 


x k (q) 


/ (2q +1 )tt \ 

x k ( q ) = J r(r)cos^w c f + — — p(t - kT) J dt 


50 


which solution is the ML estimate of the modulation. 

Making large SNR approximations, a simplified solution 55 
for the ML estimate of modulation index for M-PSK is: 


cot ft = 



r(r)sinw c r dt 


k-i 

k=0 


(7) 


60 


cos p = 


(10) 


2[(X«^(r)dr) 2 + (X^Wdr) 2 ; 


KT[(J 0 K \yM-y c «-T)?dt) + (f 0 KT (ys(r) -y s (t- T))dt) ] 


In the absence of noise this estimator produces the true 
value of modulation index. 

In yet another embodiment, MIM 110 performs modula- 
tion index estimation for the case when the modulation type 
symbol timing and data rate are also unknown. The modula- 
tion index estimators discussed above do not require an SNR 
estimate, or both an SNR estimate and a carrier phase esti- 
mate. However, they both utilize knowledge of the PSK 
modulation size, data rate, and symbol timing, as seen by the 
use (either explicitly or implicitly) of the parameters M and T 
and precise integration limits in equations (9) and (10). 

The same architecture shown in FIG. 2 above can be used 
for a general M-PSK modulation where modulation type M is 
unknown, and where the data rate (1/T) and fractional symbol 
timing (e) are also unknown. 

FCM 120, in one embodiment, includes an open frequency 
estimator when a residual carrier is available. In another 
embodiment, FCM 120 includes a closed-loop frequency 
estimator when a residual carrier is available. Notably, the 
invention contemplates use of both an additive white Gauss- 
ian noise (AWGN) channel and a Rayleigh fading channel. 

An open loop estimation of a frequency offset is obtained 
as: 


For low SNR, small argument approximations are sup- 
plied: 


A / = argmax A(A /) 
a/ 


( 11 . 1 ) 


sin h{x}=x, cos 1 


( 8 ) 
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where Af is the offset frequency and 


A(A /) = 


J] f c[k]e 


r j2xfkT s 


is the ML metric for estimating the frequency offset, and r c [k] 
is the sequence of observed baseband samples. 

However, this operation is equivalent to obtaining the fast 10 
Fourier transform (FFT) of the received sequence, taking its 
magnitude, and then finding the maximum value. 

FIG. 5 is a block diagram of an exemplary embodiment of 
open loop frequency estimator 1110. The error signal for a 
closed-loop estimator can be obtained as: 15 


e = 


dAf 


MAf) 


( 11 . 2 ) 


The derivative of L(Af) for a small € can be approximated as: 20 


e=\Y(Af+e)\-\Y(Af-e)\ 

where 


(11.3) 


25 


n- 1 (11.4) 

Y(Af + e) = 

k = 0 


An error signal detector for a closed-loop frequency cor- 
rection can be implemented based on the above equations. 

FIG. 6 is a block diagram of one exemplary embodiment of 
closed-loop frequency estimator 1220 in accordance with one 
exemplary embodiment of the invention. Notably, the dashed 
box in FIG. 6 and all other figures represent the fact that the 
hard limiter is optional. This means that closed-loop estima- 
tor 1220 may be implemented either with or without such a 
hard limiter. 

FIG. 7 is a block diagram illustrating one exemplary 
embodiment of an error signal detector 1225 for a residual 
carrier based on the above equations, wherein 


a=e J2 ™ Ts . 

FIG. 8 is a block diagram of an exemplary embodiment of 45 
an error signal detector 1225" in accordance with one exem- 
plary embodiment of the invention. Implementation of error 
signal detector 1225 M is complex. Instead of approximating 
the error signal as in (70), we note that: 


N - 1 

e = Im(r U) = 2 XN/2,N-l) 

i = 0 


( 12 ) 


where 


55 


« (13) 

X mJ , = Y j r c [k]e~> MikT ’ 

k=m 

60 

One embodiment of a closed-loop frequency estimator 
1230 for a residual carrier signal with the approximate error 
signal detector given by (12) is illustrated in the block dia- 
gram of FIG. 9. The parameters N (the number of samples to 
be summed) and 6 (gain) should be optimized and updated 65 
after the initial start to perform both the acquisition and track- 
ing of the offset frequency. 
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The gain 6 that was shown in the closed-loop frequency 
tracking system is usually part of a digital loop filter 1250. 
The digital loop filter without gain 6 can be represented as: 


One embodiment of a circuit for a digital loop filter 1250 is 
shown in the block diagram of FIG. 10. In addition to the gain 
6, the parameter b should also be optimized for achieving the 
best performance. 

Consider a data modulated signal with a suppressed carrier, 
wherein let r(t) be the received complex waveform, and a z the 
complex data representing an M-PSK or a QAM modulation. 
Let p(t) be the transmit pulse shaping. Then the received 
signal can be modeled as: 


r{t ) = ^ aiPit- iT)c J(2KAft+6c) + h{t) 


(15) 


where T is the data symbol duration, and n(t) is the additive 
white Gaussian noise (AWGN) 

with single-sided power spectral density 2N 0 W/Hz per 
dimension. 

The open loop estimation: 


Af = argmax A(A /) 
a / 


(16) 


is equivalent to multiplying the received signal by e _/(2jtA/?) , 
passing it through the matched filter with impulse response 
p(-t), and sampling the result at t=(k+ 1 )T, which produces the 
sequence of z^s. Next, sum the z^s, take its magnitude, and 
then find the maximum value by varying the frequency Af 
between -kf max and Af max where Af max is the maximum 
expected frequency offset. In one exemplary embodiment, 
FCM 120 utilizes an open-loop frequency estimator 1260 
(see the block diagram of FIG. 11) for a suppressed carrier 
with known data to perform the operations of (16). 

The error signal, for closed-loop tracking can be obtained 
as: 


An error signal detector 1275 for the closed-loop fre- 
quency correction using the above equations, in an exemplary 
embodiment, is shown in FIG. 12. Notably, in FIG. 12, DAC 
denotes a Digital -to -Analog Converter. 

The error signal can be calculated by taking the derivative 
of X 2 (Af)=IYI 2 to obtain: 


e=Im{Y*U) (18) 

The error signal in (18) can also be written as: 


e=\Y-jU\-\Y+jU\ (19) 

FIG. 13 is a block diagram of one embodiment of a closed- 
loop frequency estimator 1270 using the error signal detector 
given by (18). Similarly, a block diagram of a closed-loop 
frequency estimator 1280 using the error signal detector 
given by (19), in one embodiment, is shown in FIG. 14. 

Closed-loop frequency estimators 1270 and 1280 each 
include mixed analog and digital circuits. An all-digital ver- 
sion of a closed-loop frequency estimator 1290 in accordance 
with another exemplary embodiment of the invention is illus- 
trated FIG. 15, wherein closed-loop frequency estimator 
1290 operates on the received samples r k . In FIG. 15, p^ 
represents the discrete time version of the pulse shaping p(t), 
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where there are n samples per data symbol duration T. Nota- 
bly, an all digital version of other closed-loop estimators can 
be obtained in a similar manner and are each contemplated by 
the invention. 

FCM 120 includes a non-data-aided frequency estimator 
1292 for situations where the data is unknown to radio 
receiver 100. The I 0 (x) function can be approximated as: 




For open loop estimation, the equation: 


A/ = argmax A(A/) 
a / 


(20) 


is equivalent to multiplying the received signal by 
passing it through a matched filter with impulse response 
p(-t), and sampling the result at t=(k+l )T, which produces the 
sequence of z^s. Next, take the magnitude square of each z k , 
perform summation, and then find the maximum value by 
varying the frequency Af between - 


4 ■J'juA* 


is the maximum expected frequency offset. 

FIG. 16 is a block diagram of a circuit 1292 to perform the 
operations of (20). 

The error signal for closed-loop tracking can be approxi- 
mated as: 35 


e = 2 + ^l 2 - - ^l 2 ) 


One exemplary embodiment of an error signal detector 
1294 for closed-loop frequency correction using the above 
equations is illustrated in the block diagram of FIG. 17. 

The error signal is obtained as: 45 


e = J] Im {4(A/)w*(A/)l 


and 
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in accordance with one exemplary embodiment of the inven- 
tion. Similarly, a block diagram of one embodiment of a 
closed-loop frequency estimator 1296 using the error signal 
detector given by (23) is shown in FIG. 19. Closed-loop 
frequency estimators 1295 and 1296 include mixed analog 
and digital. 

All digital versions of other closed-loop estimators can be 
obtained in a similar manner and are contemplated by the 
invention. 

to Referring again to FIG. 1, after the modulation index is 
estimated, the incoming signal is transmitted to JEM 130. In 
accordance with one exemplary embodiment of the invention, 
JEM 130 includes a data format/pulse shape module 132 to 
determine the waveform shape and pulse of the incoming 
15 signal. In another embodiment, JEM 130 includes a signal - 
to-noise ratio (SNR) module 134 to determine the strength of 
the incoming signal relative to background noise. JEM 130, in 
yet another exemplary embodiment, includes a data rate mod- 
ule 136 to determine the data rate of the incoming signal. In 
20 still another exemplary embodiment, JEM 130 includes a 
coarse symbol synchronizer (CSS) 138 to provide symbol 
timing to various components of the radio receiver. In accor- 
dance with one exemplary embodiment of the invention, JEM 
130 includes at least two of data format/pulse shape module 
132, SNR 134, data rate 136, and symbol synchronizer 138. In 
accordance with another exemplary embodiment, JEM 130 
includes at least three of data format/pulse shape module 132, 
SNR 134, data rate 136, and symbol synchronizer 138. In yet 
another exemplary embodiment, JEM 130 includes each of 
30 data format/pulse shape module 132, SNR 134, data rate 136, 
and symbol synchronizer 138. Furthermore, each of data 
format/pulse shape module 132, SNR 134, data rate 136, and 
symbol synchronizer 138 perform their respective functions 
at substantially the same time. 

Date format module 132, in one exemplary embodiment, is 
configured to determine whether the incoming signal 
includes non-return to zero (NRZ) encoding or Manchester 
(21 A) encoding. Furthermore, date format module 132 is configured 

to determine NRZ encoding or Manchester encoding for 
40 Binary Phase-Shift Keying (BPSK) signals and/or Quadra- 
ture Phase Shift Keying (QPSK) signals, and whether these 
signals are coherently classified (i.e., one encoding format is 
always used on a specific carrier type (e.g., NRZ is always 
used on a suppressed carrier)) or non-coherently classified 
(i.e., independent of the encoding format, the modulations are 
assumed to be a fully suppressed carrier). 

In accordance with one exemplary embodiment of the 
(21B) invention, date format module 132 includes a coherent data 
format classifier 1321 for determining NRZ encoding or 
50 Manchester encoding for an incoming BPSK signal with a 
low SNR. Again, supposing the incoming signal is repre- 
sented in passband by: 


«,-(A/) 


-f 


r(t)tp(t-iT)e- j(27rAfi) dt 


( 22 ) 55 


r(f) = 2V2P ^ c n p(t-nT) cos a) c t + n(t) 


(24) 


Note that the error signal in (103) can also be written as: 


n-i (23) 

e = ^ (A/) - ju k (A/)| 2 - | Zk (A/) + ju k (A/)| 2 } 

k = o 

FIG. 18 is a block diagram of a closed-loop frequency 
estimator 1295 using the error signal detector given by (2 1 B) 


where P is the signal power, {c„} is the sequence of binary 
60 independent, identically distributed (iid) data taking on val- 
ues ±1 with equal probability, p(t) is the pulse shape (the item 
to be classified), co c is the radian carrier frequency, 1/T is the 
data (symbol) rate andn(t) is a bandpass additive white Gaus- 
sian noise (AWGN) source with single-sided power spectral 
65 density N 0 W/Hz. Then a classification choice between the 
two pulses shapes based on the LLF would be to choose 
Manchester if: 



15 
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( 2 V 2 P ) ( 2 V 2 P 

In cosh - ■ — rjfc(l) I < y In coshl r k ( 2) 


where 


(25) 


5 


* _ fti-i 

B l max~ 0— 1)AU ^ - 1 

For low SNR, (25) simplifies to: 


K - 1 

J^/c+l/2)T 


r(t) cosoj, 


M,r 

J(k+i 


r(r)cosn> c r dr < 0 


For high SNR, (25) reduces to: 


10 


15 

( 26 ) 


20 


Elf 


r(r)cosw ( 


> c tdt+ r 

J(k-\ 

£ir 


r(T)cos(x) c tdt\ < 


r(t)cos(x) c tdt - 


r(k- 

J(k+l 


( 27 ) 


r(t)cosco c t dt\ 


30 


Note that while the optimum classifier of (25) requires 
knowledge of SNR, the reduced-complexity classifier (27) 
does not. 

FIG. 20 is a block diagram of one embodiment of a low 35 
SNR classifier 1321 defined by (26) and FIG. 21 is a block 
diagram of an embodiment of a high SNR classifier 1322 
defined by (27). 

For QPSK modulation, low SNR, the classification deci- 
sion would be based on the inequality: 40 

sin [{M/2 )< |>] (28) 

FIG. 22 is a block diagram one exemplary embodiment of 
a coherent data classifier 1323 for QPSK modulation having 
low SNR. 45 

For high SNR, the following equation is equivalent to (27) 
for the non-coherent high SNR case: 


Eir' 




rtt- 

{ dt+ I 

J(k+i 

ei r 


( 30 ) 


r(t)e jujci dt\ < 


r{t)e M 


'*-r 


r{t)e W dt\ 
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FIG. 23 is a block diagram of one exemplary embodiment 
of a low SNR classifier 1324 defined by (29) and FIG. 24 is a 
block diagram of one exemplary embodiment of a high SNR 60 
classifier 1325 defined by (30). 

The solution to the problem of making an ML decision on 
the pulse shape of a modulation from a variety of different 
possibilities, in principle, follows the identical procedure as 
discussed above for data format classification, except for the 65 
fact that they are no longer restricted to digital pulse wave- 
forms. 
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Of the many measures that characterize the performance of 
a communication receiver, signal -to-noise ratio (SNR) is per- 
haps the most fundamental in that many of the other measures 
directly depend on its knowledge for their evaluation. Thus, it 
is desirable that the estimation of SNR take place with as little 
known information as possible regarding other system 
parameters such as carrier phase and frequency, order of the 
modulation, data symbol stream, data format, and the like. In 
one exemplary embodiment of the invention, SNR module 
134 is configured to estimate and/or determine the signal-to- 
noise ratio with as little known information (e.g., carrier 
phase and frequency, order of the modulation, data symbol 
stream, data format, and the like) as possible. In accordance 
with an aspect of one exemplary embodiment of the inven- 
tion, SNR module 134 includes a modified split- symbol 
moments estimator (SSME) 1341, wherein FIG. 25 is a block 
diagram illustrating an embodiment of SSME 1341. 

FIG. 26 is a block diagram of a complex baseband SSME 
1342 obtained by replacing the half-symbol accumulators in 
FIG. 25 by half-symbol I&Ds. 

Besides being used for phase compensation of the samples 
or half-symbols that enter into the expressions for computing 
U*, the frequency estimate also enters into play in determin- 
ing the estimates h* that are computed from h ± by replacing 00 
with its estimate 00 . Thus, the performance of SSME 1341 and 

1342 will depend on the accuracy of the frequency estimate 00 

with or without phase compensation. In the most general 
scenario, a taxonomy of cases for analysis is illustrated by the 
tree diagram of FIG. 27 for the sampled version of the SSME 
1341. FIG. 27 begins at the square node in the middle and 
proceeds outward to any of the eight leaf nodes representing 
combinations of oo, oo, and oo^. For SSME 1342, a few of 

the tree branches of FIG. 27, namely, 2c and 3c, do not apply. 

SSMEs 1341 and 1342 assume that the symbol timing is 
either known or could be estimated perfectly. In another 
exemplary embodiment, the SNR estimator includes SSME 

1343 (FIG. 28) to determine the SNR when symbol timing is 
imperfect but the carrier frequency is known (but the carrier 
phase is still assumed unknown). 

A generalized SSME structure that provides improved per- 
formance in the sense of lowering the variance of SNR mod- 
ule 134 will now be disclosed. To simplify matters, consider 
the ideal case of no frequency uncertainty and SSME 1342. 
Suffice it to say that the generalization is readily applied to 
SSME 1341 in a manner readily attainable by one of ordinary 
skill in the art. A generalized SSME is such that for a fixed 
observation time the quality of SNR module 134 continues to 
improve with increasing SNR. 

Even though the optimized generalized SSME requires (in 
principle) very precise prior knowledge of the true value of R, 
its performance can be reasonably well approximated by that 
of a robust estimator R^i using a very coarse prior estimate of 
R. 

Similar to the generalized SSME R L is based on complex- 
valued samples, the generalized SSME R L based on real- 
valued samples may be optimized with respect to its asymp- 
totic performance expressions 

FIG. 28 is a block diagram of one exemplary embodiment 
of time-multiplexed SSME 1343. 

The equations defining both the time-multiplexed estima- 
tor R l (SSME 1343) and the underlying observables U ± in 
terms of standard chi- squared random variables are identical 
in form to those obtained for the special case of uniform 
sub -sampling of all the symbols, L=(L, L, . . . , L). 

FIG. 29 is a flow diagram of one exemplary embodiment of 
an adaptive scheme 1344 modeled after the robust version of 
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the generalized SSME R L discussed above, wherein the inte- 
ger values of L are restricted to the set b 1= 0, 1, 2, 3, ... for 
some integer base b. 

Thus far, the behavior and performance of the SSME for 
the class of M-PSK (Mi=2) modulations with and without 
frequency uncertainty have been discussed. In various other 
exemplary embodiments, SNR module 134 includes an 
SSME structure to provide SNR estimation for offset QPSK 
(OQPSK) as well as non-constant envelope modulations such 
as QAM. As before, the performance of the estimator is still 
independent of the data symbol sequence as well as the carrier 
phase and allows for the same enhancement by increasing the 
number of pairs of observables per symbol in accordance with 
the true value of SNR. 

For the case of M-PSK, it was noted that the kth transmitted 
complex symbol in the interval (k-1) T^t^kT can be repre- 
sented in the form d k =efi k where § k takes on one of M phases 
uniformly spaced around the unit circle. A special case of the 
above corresponding to M=4 results in conventional QPSK. It 
is well-known that on nonlinear channels OQPSK provides a 
performance advantage since it reduces the maximum fluc- 
tuation in the signal amplitude by limiting the maximum 
phase change to 1 35° rather than 1 80°. Since for OQPSK, the 
complex representation of a symbol extends over 1 Vi symbols 
(because of the offset between the I and Q channels), it cannot 
conveniently be represented in the polar form d k =d^ k as 
above. Rather, the I and Q channel modulations are consid- 
ered separately. For convenience, SSME 1342 is discussed, 
however, the same conclusions also apply to SSME 1341. 

For implementing a universal SSME structure 1349 that 
will handle OQPSK as well as conventional QPSK, the results 
of inputting the I and Q channel baseband signals to half- 
symbol I&Ds are partitioned into even and odd outputs. 

Date rate module 136, in one exemplary embodiment, is 
configured to determine and/or estimate the data rate of the 
incoming signal. Knowledge of this parameter is important in 
carrying out ML detection of other parameters such as the 
carrier phase or modulation type. Here, the data rates are 
assumed to come from a set of known values and related by 
integer powers of an integer base B. This assumption allows 
estimation of the true data rate based on estimates of the 
signal-to-noise ratio (SNR) computed for various assumed 
data rates. In one embodiment, the SNR utilized is SSME 
1341 and/or SSME 1342 discussed above. These estimators 
are appealing in that the only parameter required for their 
operation is the assumed data rate. Hence, estimation of the 
data rate can be done jointly with that of the SNR. 

Though this approach provides a way of estimating both 
the data rate and SNR together, it will be shown that it is 
sensitive to symbol timing error or jitter. In fact, the presence 
of symbol timing error can severely degrade the performance 
of this estimator. To overcome this, a modification is proposed 
in which the jitter is quantized and estimated alongside the 
data rate and SNR. This approach, based on a so-called gen- 
eralized likelihood ratio test (GLRT), is robust in the presence 
of symbol timing error and can be used to jointly estimate the 
data rate, SNR, and symbol timing error all at once. The 
estimates of the symbol timing error obtained can then be 
used as coarse initial estimates for the data transition tracking 
loop (DTTL), which is used later in the receiver to obtain a 
fine estimate of the timing. 

The baseband signal received at the autonomous radio is 
assumed to consist of a constant amplitude digital data stream 
corrupted by artifacts due to the conversion from IF to base- 
band as well as additive noise. Mathematically, the received 
signal r(t) is assumed to have the following form in the com- 
plex baseband representation: 


m = J] 4 - (* + z)T)je j( “ rt+dc) + m 

5 

Here, the following is true: 

A =signal amplitude; 

d fr =kth data symbol (typically assumed to be an M-PSK 
symbol); 

10 p(t)=data pulse shape (typically either an NRZ or a 
Manchester pulse); 

T=symbol period of the data; 

€=symbol timing error (jitter) (assumed to be uniform over 
1 5 the interval [0,1)); 

coi^residual frequency offset after demodulation; 

0 c =carrier phase (assumed to be uniform over the interval 
[0,2 jt)); and 

n(t)=complex additive white Gaussian noise (AWGN) 
20 whose real and imaginary parts are uncorrelated, zero 
mean processes with power spectral density (psd) 

No 

25 2 ’ 

Prior to estimating parameters such as the carrier phase 0 C 
or the frequency offset co r , the data rate given is estimated by 
R^l/T. 

30 In one embodiment the SNR utilized is SSME 1341 for 
estimating r (t) from (32) for the case of a rectangular NRZ 
pulse shape. (Notably, for different pulse shapes, the only 
thing that needs to be changed is that the half- symbol inte- 
35 grate & dump (I&D) circuits need to be replaced with half- 
symbol matched filters) Here, T denotes the assumed symbol 
period of the system (i.e., the sample period), N denotes the 
number of system observations, and 00 , h + and h“ denote 
frequency and phase compensation factors as described 
40 above. 

The data rate estimation algorithms above are used to 
jointly estimate the data rate, SNR, and symbol timing error. 
Even with a coarse quantization of the symbol timing error, a 
rather robust estimation of the data rate in the presence of 
45 jitter is obtained. Once a reliable estimate of the data rate has 
been made, the DTTL can then be used to obtain a finer 
estimate of the symbol timing error. Furthermore, the coarse 
estimate of the jitter can be used as an initial condition for the 
DTTL which may reduce the computation time required for 
50 convergence. 

As the data rate of the received signal is not known a priori, 
at radio receiver 100 the symbol timing error corresponding 
to one specific data rate is one thing that can be independently 
55 quantized. By quantizing the jitter corresponding to one data 
rate, the quantized jitter for the remaining rates may be auto- 
matically determined. Furthermore, SSME 1341 is also a 
generalized structure to account for the quantized symbol 
timing error, and to account for computing an ensemble aver- 
60 age of the observed SNRs. An ensemble average of the 
observed SNRs is important in estimating the mean of the 
SNR of the SSME system. 

To incorporate the estimation of the quantized symbol 
timing error, JEM 130 includes CSS 138 which is described in 
65 the following algorithms. 

1 . Calculate the sequence y m from FIG. 3 0 over the range of 
values given; 
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2. Run the SSME of FIG. 26 at the highest data rate 

R s =B Imax R b . Calculate jliR^ for all n and define 
n e=argmax \i &j „ and Set i=l; 

3. Lower the assumed data rate by one step, i.e., set 


1 

Rs,new = ^s,old ■> 

and run the SSME. Calculate 


for all n and define 


4. If , increment i by 1 and go to Step 3. Otherwise, 

estimate the data rate, SNR, and symbol timing error as 
follows: 

R=B lmax ~ {i ~ l) R b 


il l - 1 

£_ - 1 

35 

As discussed above, for each assumed data rate, the SSME 
is run for each value of the quantized symbol timing error. The 
SNR and jitter for that data rate are then estimated to be the 
largest SNR and the jitter value leading to this maximum 
SNR. This data rate estimation technique halts as soon as the 40 
condition =pg. is not satisfied. This may lead to a prema- 
ture termination of the algorithm. To prevent a premature 
halting of the algorithm, a GLRT-type modification is uti- 
lized. This GLRT-type estimation algorithm is based on the 
principle that the true data rate and symbol timing error 45 
should yield the largest value of the mean of the SNR. Incor- 
rect values of these quantities, on the other hand, should lead 
to a degraded estimate of the SNR mean. As opposed to the 
previous algorithm, which lowers the assumed data rate until 
the SNR decreases, this algorithm computes the SNR for all 50 
data rates and all jitter values. The advantage to this is that it 
can prevent the algorithm from prematurely terminating, 
which can easily happen when the true SNR is low. 

With reference to data rate module 136, carrier synchroni- 
zation techniques have been developed assuming that the 55 
modulation format and signal constellation characteristics are 
known a priori. In other words, the modulation index is cho- 
sen so that either the carrier is fully suppressed or a residual 
carrier component remains. Notably, constellation character- 
istics refers to the shape of the constellation (e.g., a circle for 60 
multiple-phase shift keying (M-PSK) or a square for quadra- 
ture amplitude modulation (QAM)), and its size in terms of 
the number of signal points it contains. Aside from knowing 
the modulation index and signal constellation structure, it is 
also customary to have knowledge of the data rate and type 65 
(e.g., NRZ versus Manchester code) since the true optimum 
design of the loop depends on this information. 


20 

For radio receiver 100, the most optimistic situation would 
be that the receiver contains a carrier synchronization struc- 
ture that is capable of tracking the carrier phase independent 
of the above-mentioned considerations. While this is not 
completely possible since, for example, a squaring loop (or 
equivalently a BPSK Costas loop) cannot track a QPSK 
modulation and likewise a 4 th power loop (or equivalently, a 
QPSK Costas loop, sometimes referred to as an inphase- 
quadrature (I-Q) loop) cannot properly track a BPSK signal. 
Nevertheless, while in principle each carrier synchronization 
loop developed for a given modulation format, constellation, 
and data rate/type has certain unique characteristics, they do 
share a number of similarities (e.g., a common front end 
demodulator structure, that allows one to consider designs 
that could be operational in the absence of complete a priori 
knowledge of all of these characteristics). For example, if the 
modulation is restricted to the M-PSK class, then it is possible 
to construct a universal structure that performs the carrier 
synchronization function for all values of M. This structure is 
derived by first determining the maximum a posteriori (MAP) 
estimate of carrier phase based on an observation of the 
received signal, namely, M-PSK plus additive white Gaussian 
noise (AWGN) and then using this to motivate a closed loop 
carrier synchronization loop. By making an analogy between 
the closed loop bandwidth and the noise bandwidth of an 
integrate-and-dump (I&D) filter of duration equal to the 
observation time for the open loop MAP estimate, the closed 
loop approaches the Cramer-Rae lower bound on the variance 
of an unbiased estimate of the phase of a modulated carrier. 

Still further, if the modulation is known to be other than 
suppressed carrier (i.e., a modulation index less than jt/2 
radians), then it is still possible to exploit the power in both 
the data and residual carrier components for carrier tracking 
purposes provided one has knowledge of the modulation 
index itself. Such knowledge could be derived non-coher- 
ently (i.e., in the absence of carrier synchronization) from a 
suitable modulation index estimator. Loops of this type are 
often referred to as “hybrid carrier tracking loops,” and like 
their suppressed carrier counterparts, are motivated by the 
same MAP considerations. 

In various exemplary embodiments, the matched filters in 
the I and Q arms of the loop are replaced with simple low pass 
filters to make the carrier synchronizer operation somewhat 
independent of the exact pulse shape, but without some atten- 
dant loss in performance. These systems include a fixed 
modulation bandwidth which implies a fixed data symbol rate 
for all values of M. These systems include simplifications 
based on low SNR approximations applied to the nonlineari- 
ties inherent in the MAP phase estimate. When this is done, 
the error signal in the loop for M-PSK is of the form sin M(|), 
where <|) is the loop phase error, which from simple trigonom- 
etry can be written as sin M$=2 sin [(M/2)(|)] cos [(M/2)<|>]. 
Thus, it is seen that the error signal in the loop for M-PSK is 
formed from the product of the error signal sin [(M/2)(|)] and 
the lock detector signal cos [(M/2)(|)] in the loop for M/2-PSK 
modulation. This relationship forms the basis for implement- 
ing a universal structure. 

In the past, carrier synchronization loops have typically 
fallen into two categories: those that track a discrete carrier 
(e.g., the phase-locked loop (PLL)) and those that track a fully 
suppressed carrier (e.g., the Costas loop). A fully suppressed 
carrier comes about when a digital modulation is impressed 
on a carrier with a modulation index equal to jt/2 radians, 
whereas a discrete (residual) carrier component appears in the 
spectrum when the modulation is less than jc/2 radians. 

For implementing a universal SSME structure that will 
handle OQPSK as well as conventional QPSK, the results of 
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inputting the I and Q channel baseband signals to halfsymbol 
I&Ds are partitioned into even and odd inputs. SNR estimator 
134 can use the other embodiments SSME 1341-1344 for 
estimating SNR when QAM is transmitted. 

Date rate module 136, in one exemplary embodiment, is 5 
configured to determine and/or estimate the data rate of the 
incoming signal. Knowledge of this parameter is important in 
carrying out ML detection of other parameters such as the 
carrier phase or modulation type. Here, the data rates are 
assumed to come from a set of known values and related by 
integer powers of an integer base B. This assumption allows 
estimation of the true data rate based on estimates of the 
signal-to-noise ratio (SNR) computed for various assumed 
data rates. In one embodiment, the SNR utilized is SSME 15 
1341 and/or SSME 1342 discussed above. These estimators 
are appealing in that the only parameter required for their 
operation is the assumed data rate. Hence, estimation of the 
data rate can be done jointly with that of the SNR. 

Though this approach provides a way of estimating both 20 
the data rate and SNR together, it will be shown that it is 
sensitive to symbol timing error or jitter. In fact, the presence 
of symbol timing error can severely degrade the performance 
of this estimator. To overcome this, a modification is proposed 
in which the jitter is quantized and estimated alongside the 25 
data rate and SNR. This approach, based on a so-called gen- 
eralized likelihood ratio test (GLRT), is robust in the presence 
of symbol timing error and can be used to jointly estimate the 
data rate, SNR, and symbol timing error all at once. The 
estimates of the symbol timing error obtained can then be 30 
used as coarse initial estimates for the data transition tracking 
loop (DTTL), which is used later in the receiver to obtain a 
fine estimate of the timing. 

To account for the presence of symbol timing error, typi- 
cally a digital transition tracking loop (DTTL) is used. How- 35 
ever, the DTTL architecture requires knowledge of the data 
rate to function as required. To overcome this restriction, the 
proposed SSME based data rate estimator can be modified to 
provide a coarse estimate of the symbol timing error in addi- 
tion to the data rate, which can then be used to obtain a finer 40 
estimate of the symbol timing error using a DTTL. As the data 
rate of the received signal is not known a priori, at radio 
receiver 100 the symbol timing error corresponding to one 
specific data rate is one thing that can be independently quan- 
tized. By quantizing the jitter corresponding to one data rate, 45 
the quantized jitter for the remaining rates may be automati- 
cally determined. 

As the number of quantization steps increases exponen- 
tially as the assumed data rate decreases, it is tempting to 
think that a better estimate of the data rate, SNR, and symbol 50 
timing error for lower true data rates than for higher rates will 
be obtained. However, this is offset by the fact that for a fixed 
observation time interval, an exponentially larger number of 
observations for higher true data rates than for lower ones is 
obtained. Hence, there is an implicit tradeoff between the 55 
number of signal observations and the number of jitter quan- 
tization levels for each true data rate. 

One of the advantages of uniformly quantizing the symbol 
timing error to N e6ySr steps is that it leads to an efficient all 
digital implementation of the SSME system. 60 

Furthermore, SSME 1341 is also a generalized structure to 
account for the quantized symbol timing error, and to account 
for computing an ensemble average of the observed SNRs . An 
ensemble average of the observed SNRs is important in esti- 
mating the mean of the SNR of the SSME system. To do this, 65 
the discrete-time signal y m is partitioned into blocks over 
which the SNR is to be computed. For each block, the SSME 


computes an estimate of the SNR and then an ensemble 
average of the SNR is computed over the blocks. 

There are several things to note regarding the structure of 
all digital SSME system 1362. First, all digital SSME system 
1362 includes discrete summations analogous to the sampled 
version of the SNR estimator discussed above with regard to 
SNR 134. Furthermore, note that all of the signals starting 
from the half signal integrator outputs are indexed with a 
semicolon followed by n. Finally, note that to form a single 
SNR estimate, a total of B^N^, samples are squared and 
accumulated. This was chosen as such here to keep the total 
observation time interval or epoch block fixed. 

All other things being equal, a modulation classifier that 
requires the least knowledge of channel parameters is pre- 
ferred. If a classifier required only as, for example, it would be 
preferred over one requiring P, and a 2 separately. 

FIG. 32 is a block diagram of one exemplary embodiment 
of a modulation classifier 1410. It is the qLLR modulation 
classifier. Beginning with a quasi-log-likelihood ratio 
(qLLR) approximation, and the further approximation 
I 0 (x)se x , valid for large x, (high post-detection SNR) leads to 
the approximate threshold: 


N(2P t ) M/2 ( 34 ) 


This approximate threshold actually outperforms the 
“optimum” threshold for symbol SNR’s above 4 dB. 

Alternatively, in another exemplary embodiment the qLLR 
metric 1410 may be modified by normalizing its magnitude, 
resulting in the nqLLR metric: 


N - 1 

nqLLR = — 

N ~ l M 
£ \fn\ M 
n= 0 

The analytical derivation of appropriate thresholds for the 
qLLR in (34) depended on the particular form of the metric, 
and involved some approximation. As an alternative, an 
empirical method to optimize the threshold can be used for 
any classifier metric, including the qLLR and nqLLR metrics 
is disclosed. Empirical threshold determination consists of 
the following process: 

1 . Generate a large number of received noisy M-PSK and 
M'-PSK samples; 

2. Group the samples into blocks of length N, and compute 
a set of sample classifier metrics for the qLLR classi- 
fier); 

3. Sort the metrics for M-PSK and M'-PSK received sym- 
bols separately, in order of increasing value; 

4. For each observed metric of value a, in increasing order: 

(a) Count the number of observed M-PSK metrics hav- 
ing value less than a, 

(b) Count the number of observed M'-PSK metrics hav- 
ing value greater than a, and 

(c) Compute the probability of misclassification when 
using threshold a, using the sum of the above counts; 
and 

5. Report the threshold that minimizes the probability of 
misclassification. 

This procedure may be implemented efficiently enough 
that a desktop machine can generate about 100 million 
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samples and determine an empirically optimum threshold in 
less than an hour. The following additional points should be 
noted: 

1 . Threshold optimization need be performed once, offline. 

In a practical implementation, given N and a table look-up 5 
may be used to determine the threshold. 

2. The optimum threshold for the LR metric and its 
approximation in (344) is zero, and requires no empirical 
optimization. 

3. The optimum threshold for the qLLR metric for BPSK/ 10 
QPSK classification is relatively flat over a broad region of y s 
being near 0.6 for y s >-5 dB. This is a desirable characteristic, 
because may not be known exactly. 

4. The optimum threshold for the nqLLR metric is also 
quite flat, and ranges from about 0. 1 5 to 0.6 for - 1 0<y 5 .<l 0. At 15 
high SNR, it is approximately equal to the optimum qLLR/ 
(2NPj threshold, because the qLLR/(2NP ? ) metric and the 
nqLLR metric are nearly identical in that region. Thus, the 
nqLLR and qLLR classifiers have similar performance, and 
the nqLLR has the advantage that neither P ? nor a 2 need be 20 
known. 

Thus far exact and approximate likelihood-based non-co- 
herent classifiers of M-PSK signals have been discussed 
under the assumption of unknown (uniformly distributed) 
carrier phase, but with perfectly known fractional symbol 25 
timing. In other embodiments these classifiers are extended to 
the case where the symbol timing is unknown and also uni- 
formly distributed 

To simplify the implementation of the ML classifiers, in 
one embodiment, approximations of the nonlinearities 30 
involved in their definitions are used in much the same way as 
was done for the conventional (non-offset) modulations. The 
approximate LLF is: 


^ = X^Z [f/„l 2 + lfe„l 2 ] + ln/o TZa\Yj ti -tQj 


35 


(35) 


For BPSK), the following LLF is obtained: 40 


tan hx^sgiLX 


(36) 


where again x„ (0; 0 c )=x /w (0j since in this case x Qn (0 C )=O. 

To get to the final simplification, the approximation \ n 
I 0 (x)=x (x»1) is applied to result in the qLLR: 45 


100 employs symbol timing estimation and tracking schemes 
whose implementations can adapt themselves to this chang- 
ing environment using the knowledge obtained from SNR 
module 134. Furthermore, since SNR module 134 itself 
requires knowledge of symbol timing, circuitry for obtaining 
a coarse estimate of this timing is important. 

In one embodiment, FSS 1 60 includes a DTTL is motivated 
by an open loop estimate to determine the MAP estimate of 
the symbol timing based on an observation of a block of N 
symbols. Analogous to the ML approach taken in modulation 
classification, the likelihood function (LF) of the received 
signal vector is formed conditioned on the unknown param- 
eters to be estimated 

For the purpose of finding the MAP estimate of symbol 
synchronization alone, perfect knowledge of the carrier phase 
is assumed, and wherein 0 C =O. Under this assumption, the 
MAP estimate of symbol timing £ MAP MAP is given by: 


e MAP = argmax exp 



cosh[x„(#; e)] 


(44) 


where now: 




r l+ % P (, -nT- eT)dt\ 
No { J(n+e)T J 


(45) 


Note that the actual fractional symbol timing offset € is 
embedded in the received complex baseband signal r(t) and 
thus the difference between € MAP and € represents the nor- 
malized symbol timing error. 

A closed-loop configuration that implements the expres- 
sion in (446) as an error signal is referred to as a MAP 
Estimation Loop. 

For binary signals (M=2), the error signal of (45) simplifies 
to: 


qLLR' = 


4 cr 4 


Z DM 1 - l ? eJ 2 ] + 2 Z d- Z ti. ~ 


(37) 


e = 2_j tanh [ a^(0; s)] — s ) 


(46) 


50 


The operation and performance of various receiver func- 
tions can be quite sensitive to knowledge of the timing (data 
transition epochs) of the received data symbols. Thus, the 
ability to accurately estimate this parameter and continuously 55 
update the estimate (i.e., perform symbol synchronization) 
with as little knowledge of other parameters is important to 
successful operation of a radio receiver 100 . 

Traditionally, symbol synchronization techniques have 
been developed assuming that the data symbols are binary, the 60 
modulation format (e.g., NRZ or Manchester data) is known 
a priori, and carrier synchronization is perfect. Thus, the 
symbol synchronization problem has been solved entirely at 
baseband assuming perfect knowledge of the carrier phase 
and frequency. 65 

Since radio receiver 100 functions, in general, over a wide 
range of SNRs, in one exemplary embodiment, radio receiver 


where: 


2 yfp rin+l+s) 

x n (0;s)= —— I r{t)p{t -nT - sT) d t 

No J(n+£) 


4(0; s ) = ■ 


2 VP r(n+i+£) dp{t-nT- sT) 


No 


■f 

J(n-\ 


r(t) 


dt 


dt 


and r(t) is now a real signal. A block diagram of one 
embodiment of a MAP estimation loop 1610 that uses e of 
(46) as an error signal to control a timing pulse generator is 
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illustrated in FIG. 34, wherein the shorthand notation T„(€) 
represents the time interval (n+e)T^t^(n+ 1 +€)T. In FIG. 34, 
the accumulator represents the summation over N in (46). 
Thus, based on the above model, the loop would update itself 
in blocks of N symbols. In one embodiment, this block-by- 5 
block accumulator is replaced by a digital filter which updates 
the loop every T seconds, and whose impulse response is 
chosen to provide a desired dynamic response for the loop. 
The design of this filter and its associated closed-loop 
response characteristic are not dictated by the MAP estima- 10 
tion theory which explains the use of the term “MAP-moti- 
vated” when describing the MAP estimation loop. 

To go from the MAP estimation loop to the conventional 
DTTL, one needs to approximate the hyperbolic tangent non- 
linearity for large values of its argument, equivalently, at high 1 5 
SNR and 2) characterize, i.e., approximate, the derivative of 
the pulse shape required in (50). Specifically, for large values 
of its argument, one has the approximation: 

tan hx=sgnx (51) 

Based on the above assumptions and discussion, it is now 
clear that the MAP estimation loop migrates to a DTTL 1620 
as illustrated in the block diagram FIG. 35. 

In one exemplary embodiment, FSS 160 includes a DTTL 
that is motivated by the MAP estimate of symbol timing at 25 
low SNR and, in particular, the I arm hard decisions are 
replaced by soft decisions whereupon in the limiting case, the 
hard limiter is replaced by a linear device. Under high SNR 
conditions where the nonlinearity is approximated as in (51), 
in one embodiment, the I arm of the resulting symbol syn- 30 
chronizer becomes a detector of a transition in hard decisions 
made on successive symbols. In accordance with one aspect 
of one embodiment, the DTTL includes a synchronizer that 
results from approximating the nonlinearity for small values 
of its arguments as appropriate at low SNR. 35 

For sufficiently small values of its argument, the hyper- 
bolic tangent nonlinearity can be approximated by: 

tan hx ** sgnx (52) 

(i.e., a linear function). Under this assumption, a MAP- 40 
motivated closed-loop synchronizer 1630 is illustrated in the 
block diagram of FIG. 35, and its performance is analyzed as 
follows. After known carrier phase demodulation by the car- 
rier reference V2cos(co c t+0 C ), the baseband signal input to the 
LDTTL is given by: 45 


= s(t, e) + n{t ), s{t, s) = \j~P ^ d n p{t - nT - sT ) 


where, consistent with the assumption of NRZ data, p(t) is 
a unit amplitude rectangular pulse of duration T sec and {d„} 
is an iid±l sequence with d w representing the polarity of the 
nth data symbol. The additive noise is a white Gaussian 55 
process with single- sided power spectral density N 0 W/Hz. 
The local clock produces a timing reference for the I and Q 
I&D filters that depends on its estimate € of e. 

One can immediately apply the same small and large argu- 
ment approximations to the hyperbolic tangent nonlinearity 60 
and, analogous to FIGS. 35 and 36, arrive at DTTL-like 
implementations which are illustrated in the block diagrams 
of FIGS. 37 and 38. 

As discussed above, open loop MAP estimation of symbol 
epoch involves finding the conditional (on the symbol timing) 65 
likelihood function of the received signal based on a single 
observation of the received signal over a block of symbols. 


Furthermore, since the unknown symbol epoch is assumed to 
be uniformly distributed over the symbol interval, the MAP 
estimate is equivalent to the ML estimate. Moreover, the 
traditional closed-loop estimation scheme motivated by the 
MAP estimation approach employs an error signal derived 
from the derivative of the CLF that can be updated at intervals 
corresponding to the symbol time. Since for rectangular 
pulses (e.g., an NRZ data stream) the derivative of the CLF, 
which is related to the derivative of the pulse shape, is unde- 
fined, closed-loop structures motivated by the MAP estima- 
tion approach generally do not exist. Nevertheless, with suit- 
able approximations of the derivative of the pulse shape, such 
a closed loop (e.g., the DTTL) will in fact provide symbol 
synchronization for an NRZ data stream with rectangular 
pulses. The closed-loop approach provides a continuous 
updating (tracking) of the symbol timing (once per bit inter- 
val) which is desirable in the presence of channel dynamics 
whereas the open approach is usually regarded as either a one 
shot estimator (i.e., compute the MAP or MMS estimate 
based on a single observed long block of data), or a block- 
by-block estimator where the single shot is sequentially 
repeated over and over. 

What is important to observe is that the open loop estima- 
tion techniques can be modified to provide sequential updates 
at the symbol rate to the symbol timing epoch estimates and as 
such resemble the closed-loop techniques with however, 
improved performance. As such, the presented invention 
includes a sequential digital implementation of the MAP 
estimation of symbol epoch that can track the dynamics in this 
parameter yet provide a performance approaching that of the 
true optimum MAP estimation technique. 

Because of the assumption that the unknown parameter 
being estimated (i.e., symbol epoch) is constant over the 
observation, the implementation in FIG. 36 (referenced as 
1640) includes a running accumulator with uniform weight- 
ing. In other embodiments where the parameter is dynamic 
but slowly varying, one can assume that the unknown param- 
eter is constant over a finite number of symbol intervals. As 
such, the uniform running accumulator is replaced with a 
weighted running accumulator which reflects a fading 
memory and is analogous to what is done in a closed-loop 
architecture by using a digital filter following the error signal. 
This embodiment (referenced as 1650) is illustrated in the 
block diagram of FIG. 37, wherein the feedback term is 
multiplied by a constant a<l . This achieves a running accu- 
mulator with a geometric weighting which has the input- 
output characteristic: 




Since the running accumulator also accomplishes the data 
detection (matched filter) function, then the epoch estimate 
index, i^ of (55) may be used to determine in each symbol 
interval which running accumulator output to use for making 
a hard decision on that symbol. 

When the pulse shape is other than rectangular, the modi- 
fication that allows the bank of N c “ correlations to be replaced 
by a sliding accumulator as in FIG. 36 is not possible. How- 
ever, if N g is large and the pulse shape is approximated by a 
piecewise constant staircase function with N G steps, then the 
correlation of the received signal and the pulse shape in a 
quantization interval can be written as: 


$J + ^rii)p{i)dt=pl ,^ A r{t)dt 
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where p z is the assumed constant value of p(t) in the interval 
iA^t^(i+l)A. In view of (56), one modification of FIG. 34 
that is necessary to allow for the inclusion of an arbitrary 
pulse shape is to replace the uniform sliding accumulator with 
a weighted sliding accumulator (see FIG. 36) where the 
weights are equal to the piecewise constant values of p(t). 
Furthermore, for sufficiently large N e , one can approximately 
replace the integral of over the quantization interval by the 
value of r(t) at the midpoint of this interval times the duration 
of this interval, A. As such, the integrator at the input of FIG. 
34 can be replaced simply by a uniform sampler at rate 1/A. 

In addition to “optimum” symbol synchronizers, such as 
the ones discussed thus far that are motivated by the MAP 
estimation approach, the present invention encompasses 
other embodiments for determining the symbol timing. A 
block diagram of a real non-coherent version of a cross- 
spectrum synchronizer 1660 is illustrated in the block dia- 
gram of FIG. 38, wherein the input is a bandpass signal whose 
signal component is given by s 2 (t,€)=Re {s(t,€)e /Wc ^}, with co c 
denoting the carrier frequency. 

Another exemplary embodiment of FSS 160 includes non- 
coherent versions of a DTTL. The low SNR version of the 
noncoherent DTTL, herein given the acronym NC-LDTTL is 
a parallel combination of two independent coherent LDTTLs 
acting on the I and Q baseband signals. FIG. 39 is a block 
diagram illustrating NC-LDTTL 1670. In one exemplary 
embodiment, it is important that a coarse estimate of symbol 
timing be provided to SNR module 134, and is derived from 
the same statistics as used to form this estimator itself. 

Benefits, other advantages, and solutions to problems have 
been described herein with regard to specific embodiments. 
However, the benefits, advantages, solutions to problems, and 
any element(s) that may cause any benefit, advantage, or 
solution to occur or become more pronounced are not to be 
construed as critical, required, or essential features or ele- 
ments of the invention. All structural, chemical, and func- 
tional equivalents to the elements of the above-described 
exemplary embodiments that are known to those of ordinary 
skill in the art are expressly incorporated herein by reference 
and are intended to be encompassed by the invention. 

The invention claimed is: 

1. A method for autonomously receiving and decoding 
incoming signals of a variety of signal types without a priori 
knowledge of the defining characteristics of the incoming 
signals, comprising the steps of: 

receiving a signal of unknown signal type; 

making a first maximum likelihood estimate of a defining 
characteristic of said unknown signal type; 

utilizing said first estimate to make a second maximum 
likelihood estimate of another defining characteristic of 
said unknown signal type; 

selecting an appropriate decoder based on said first and 
second estimates; and 

decoding said signal with the selected decoder. 

2. The method for autonomously receiving and decoding 
incoming signals of a variety of signal types without a priori 
knowledge of the defining characteristics of the incoming 
signals according to claim 1, wherein said step of making the 
first maximum likelihood estimate comprises a substep of 
making a maximum likelihood estimate of a modulation 
index of said received a signal of unknown signal type. 

3. The method for autonomously receiving and decoding 
incoming signals of a variety of signal types without a priori 
knowledge of the defining characteristics of the incoming 
signals according to claim 2, wherein said substep of making 
a first maximum likelihood estimate fiirther comprises a sub- 
step of tracking an offset frequency of said received signal of 
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unknown signal type and outputting a signal having an esti- 
mated corrected frequency based on the incoming signal and 
the estimated likelihood estimate of the modulation index of 
said received a signal of unknown signal type. 

5 4. The method for autonomously receiving and decoding 

incoming signals of a variety of signal types without a priori 
knowledge of the defining characteristics of the incoming 
signals according to claim 3, wherein said step of making a 
first maximum likelihood estimate further comprises a sub- 

10 step of mixing the received signal of unknown signal type 
with the outputted signal having an estimated corrected fre- 
quency. 

5. The method for autonomously receiving and decoding 

15 incoming signals of a variety of signal types without a priori 

knowledge of the defining characteristics of the incoming 
signals according to claim 4, wherein said step of making a 
first maximum likelihood estimate further comprises a sub- 
step of deriving approximations of a maximum likelihood 

20 modulation classification based on said mixed signal from 
said mixing sub step. 

6. The method for autonomously receiving and decoding 
incoming signals of a variety of signal types without a priori 
knowledge of the defining characteristics of the incoming 

25 signals according to claim 3, fiirther comprising a step or 
tracking a carrier phase of said received signal of unknown 
signal type. 

7. A system for receiving and decoding incoming signals 
without a priori knowledge of any defining characteristics of 

30 the incoming signals types, comprising: 

a radio receiver adapted to receive an incoming signal, said 
radio receiver including, 

a first estimator for making a first maximum likelihood 
estimate of a defining characteristic of said incoming 

35 signal, 

a second estimator for utilizing the first maximum likeli- 
hood estimate to make a second maximum likelihood 
estimate of another defining characteristic of said 
incoming signal, 

40 a decoder selector for selecting an appropriate decoder 
scheme based on said first and second maximum likeli- 
hood estimates, and 

at least one decoder for implementing the selected decoder 
scheme to thereby decode said signal. 

45 8 . The system for receiving and decoding incoming signal s 

according to claim 7, wherein said receiver further comprises 
a modulation index module for making a first maximum 
likelihood estimate of a modulation index of said incoming 
signal. 

50 9 . The system for receiving and decoding incoming signal s 

according to claim 8, wherein said receiver further comprises 
a frequency correction module in communication with said 
modulation index module for tracking an offset frequency of 
said incoming signal and for outputting a signal having an 

55 estimated corrected frequency based on the incoming signal 
and the estimated modulation index from said modulation 
index module. 

10. The system for receiving and decoding incoming sig- 
nals according to claim 9, wherein said receiver fiirther com- 

60 prises a multiplier for mixing the incoming signal with the 
outputted signal having an estimated corrected frequency 
from the frequency correction module. 

11. The system for receiving and decoding incoming sig- 
nals according to claim 10, wherein said receiver further 

65 comprises a joint data rate estimator, signal -to-noise ratio 
estimator, pulse shape classifier, and symbol timing estima- 
tor. 
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12. The system for receiving and decoding incoming sig- 
nals according to claim 11, wherein said receiver further 
comprises a modulation classifier module for deriving 
approximations of a maximum likelihood modulation classi- 
fication based on said mixed signal from said multiplier and 
the outputted signal with estimated corrected frequency. 

13. The system for receiving and decoding incoming sig- 
nals according to claim 12, wherein said receiver further 
comprises a carrier synchronization module for tracking a 
carrier phase of said incoming signal. 
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14. The system for receiving and decoding any of a variety 
of incoming signals of various signal types according to claim 
13, wherein said at least one decoder further comprises a 
plurality of decoders, said receiver selecting an appropriate 
decoder based on said first and second maximum likelihood 
estimates to thereby decode said signal. 



